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Abstract 

113 residual product nuclide yields in a 1.0 GeV proton-irradiated thin monoisotopic 208p]^ 
sample and 107 residual product nuclide yields in a 2.6 GeV proton-irradiated ""*W sample 
have been measured and simulated by 8 different codes. The irradiations were made using 
proton beams extracted from the ITEP synchrotron. The nuclide yields were 7-spectrometered 
directly using a high-resolution Ge-detector. The 7-spectra were processed by the GENIE-2000 
code. The ITEP-developed SIGMA code was used together with the PCNUDAT nuclear decay 
database to identify the 7-lines and to determine the cross-sections. The ^^Al(p,x)^^Na reaction 
was used to monitor the proton flux. The measured yields are compared with calculation the 
LAHET (with ISABEL and Bertini options), CEM95, CEM2k, CASCADE, CASCADE/INPE, 
INUCL, and YIELDX codes. Estimates of the mean deviation factor are used to demonstrate 
the predictive power of the codes. The results obtained may be of interest in studying the 
parameters of the Pb and W target modules of the hybrid Accelerator-Driven System (ADS) 
facilities. 
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Introduction 

At present, the Pb-Bi eutectic and W are regarded as the most promising target materials for 
ADS facihties |j|]-|3|. As a result, the high-energy irradiation mode of using the materials ne- 
cessitates additional studies of the nuclear-physics characteristics of Pb, Bi, and W, particularly 
the yields of residual product nuclei under proton irradiation in a broad range of energies from 
a few MeV to 2-3 GeV. Results of such studies are extremely important when designing even 
demonstration versions of the ADS facilities. 

Undoubtedly, computational methods will play an important role when forming a set of nu- 
clear constants for ADS facilities. Therefore, verification of the most extensively used simulation 
codes has proved to be of a high priority. 



Basic definitions and computational relations 

The formalism of representing the reaction product yields (cross sections) in high-energy proton- 
irradiated thin targets is described in sufficient detail in [Q. In terms of the formalism, the 

variations in the concentration of any two chain nuclides produced in an irradiated target (A^i — 

N2 may be presented to be a set of differential equations that describe the production and 
decays of the nuclides. By introducing a formal representation of the time functions, Fi, of the 

_ 1 _ -XiKT 

form Fi = (1 — e~ '"^j — -X T ' °^ another monitor product; r is the 

duration of a single proton pulse; T is the pulse repetition period; K is the number of pulses 
within the irradiation period), which characterize the nuclide decays within the irradiation time, 
and by expressing (similar to the relative measurements) the proton fluence size via the monitor 
reaction cross section, agt, we can present the cumulative and independent yields as: 

Aq NAlFNa 

CTst (1) 



■qiElFiNNa Nt Xnq 



^cum ^ ^1 NaI A2 - Ai F^a ^ 

^ l'iri2£2FlNNa Nt X2 >^Na 

ind_ M2 AA 1 NAlFNa 

' V^2 F1X2) me2NNaNT XNa 

cum And. ^, cum / ^1 , ^2 \ 1 NaI F^g 

where af^^ is the cumulative cross section of the first nuclide; o"!"'^ and cr™™' are the independent 
and cumulative cross sections of the second nuclide; Nai and Nj^ are the numbers of nuclei in 
the monitor (standard) and in experimental sample, respectively; rji and r]2 are the 7-line yields; 
El and 82 are the spectrometer efficiencies at energies E^^ and E^ji ^1 is the branching ratio of 
the first nuclide; Ai, A2, and X]\fa are, respectively, the decay constants of the first and second 
nuclides and of the monitor product (■^^Na and/or ■^^Na). 

The factors Aq, Ai, and A2 are calculated through fitting the measured counting rates in 
the total absorption peaks, which correspond to energies E^^ (the first nuclide) and E^j (the 
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second nuclide), by exponential functions. It should be noted that formulas (|T|)-(§) were derived 
on assumption that the 7-intensities of the two nuclides produced under irradiation are recorded 
up to the desired accuracy within an interval of time from the irradiation end to the moment of 
the ultimate detectable intensity. If, for some reasons, the factor Ai cannot be found, then the 
factor A2 will be used together with expression (14) from |^ to determine the quantity cjg"™*, 
which we called the supra cumulative yield: 

cum* „ I "^2 cum 

CF2 = (72 + T T-Z^l<7i = 

Al — A2 

A2 Nai FNa 

-CTst (5) 



The resultant value of o"™'"* may prove to be very different from Nevertheless, the 

supra cumulative yield can be either used directly to verify the codes, or determined further 
up to (T2"™' if the needed data are obtained elsewhere (for example, from inverse-kinematics 
experiments). 



Experimental techniques 

A 10.5 mm diameter, 139.4 mg/cm^ monoisotopic ^osp^^ metal foil sample (97.2% ^'^^Pb, 
1.93% 207pi3^ 0.87% 206pb, <0.01% 204pi3^ < 0.00105% of chemical impurities) and a 38.1 
mg/cm^ ^"-^W metal foil sample (99.95% W, <0.05% of chemical impurities), both of 10.5- 
cm diameter, were proton- irradiated. 139.6 mg/cm^ and 139.1 mg/cm'^ Al foils of the same 
diameter were used as monitors. Chemical impurities of the monitor did not exceed 0.001%. 

The samples were irradiated by the external proton beam from the ITEP U-10 synchrotron 
H . The average flux densities during irradiation of Pb and W samples were of 1.4x10^° p/cm^ 
and 2.8x10^'^ p/cm^, respectively. 

Our measurements were supported by extra researches aimed at reducing the systematic 
errors in the experimental results. Theses researches included: 

• experiments to specify the neutron component in the extracted proton beams, 

• experiments to specify the ^'^Al(p,x)^^Na monitor reaction cross section, 

• studies to specify the dependence of the 7-spectrometer detection efficiency on the position 
geometry of irradiated sample, 

• studies to optimize the 7-spectrum simulation codes. 

Figs. H and |2| show the results of measuring the neutron component in the extracted 
proton beams, i.e., the neutron-to-proton flux density ratio, ^n/^p- Fig. ^presents the monitor 
reaction cross sections measured here and in other worksQ. The height-energy dependence of the 
detection efficiency is displayed in Fig. ^ 

^ MI85 - Nucl. Phys. A 441 (1985) 617; MI86 - NIM B 16 (1986) 61; MI89 - Analyst 114 (1989) 287; Mi90 
- NEANDC(E)-312-U(1990) 46; MI93 - INDC(GER)-037/LN(1993) 49; MI95 - NIM B 103 (1995) 183; MI96 - 
NIM B 114 (1996) 91; MI97 - NIM B 129 (1997) 153. 
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Q_ 1 a T O - 0.2 GeV (SSNTD) 
iQ, l\\ □ - 0.2 GeV (Ge-deteotor) ; 
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■q 1 



O - 0.8 GeV (SSNTD) 



H 



O - 2.6 GeV (SSNTD) 
i □ - 2.6 GeV (Ge-detector) ; 
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Figure 1: The neutron backgrounds around the extracted pro- Figure 2: Neutron compo- 
ton beams that irradiate thin experimental samples. nent in the extracted proton 

beams of different energies. 
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Proton energy [MeV] [keV] 



Figure 3: The ^ Al(p,x)^'^Na monitor re- 
action cross sections measured in this and 
other works. 



Figure 4: The experimental and calculated de- 
tection efficiency of the spectrometer. 
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The discrepancies between the two sets of high-energy data in Fig. ^ have yet to be studied. 
The analytical expression of the spectrometer detection efficiency as a function of energy and 
sample position height is 



eiE,H)=ebase (E) 



jqi +q2 ■ In E + Hbase) 

{qi + q2-lnE + H) 



(6) 



where qi and q2 are parameters defined by fitting the experimental results. An analysis of the 
7-spectrum processing codes has shown that the GENIE-2000 code is superior to the others 
because of its interactive mode of fitting the peaks, which permits correction of the automated 
computer-aided processing; therefore, we chose it for our work. 



Experimental results and measurement errors 

Tables |l| and ^ present the results of measuring the reaction product yields in the 1 GeV proton- 
irradiated 208p]-, 2.6 GeV proton-irradiated ""*W samples. 113 yields from ^osp]^ have been 
obtained, of which, 6 independent yields (i), 17 independent yields of metastable states (m), 
15 independent yields of metastable and ground states (Smj-|-g), 64 cumulative yields (c), and 
11 supra cumulative yields, when the addend may exceed the determination error (c*). 107 
yields from "■"*W are presented, of which, 6 independent yields (i), 9 metastable state yields 
(m), 5 yields of metastable and ground states (Smj+g), 86 cumulative yields (c), and 1 supra 
cumulative yield (c*). 

Prom Tables |l] and Q one can see that the experimental errors are ranging within ~(6-35)%. 
The main contribution to the total error is from uncertainties in the nuclear data, namely, in 
the absolute quantum yields and cross sections of the monitor reactions. 



Comparison with experimental data obtained elsewhere 

Table | and Fig. | compare some of the present results with experimental data of other 
laboratories published in [Q. 

Table 1: Experimental product nuclide yields in 1 GeV proton-irradiated ^osp^-, 



Product 


Tl/2 


Type 


Yield (mb) 




6.243d 


i 


4.60± 0.29 


205Bi 


15.31d 


i 


6.20± 0.40 




11.22h 


i(ml+m2+g) 


5.29± 0.80 


203Bi 


11.76h 


i(m+g) 


4.84± 0.59 


204mp|^ 


67.2m 


i(m) 


11. 0± 1.0 


203pb 


51.873h 


c 


31. 5± 2.1 


201pb 


9.33h 


c* 


26. 9± 2.4 


200pb 


21. 5h 


c 


18. 2± 1.2 


198pb 


2.4h 


c 


8.9± 2.1 


197mp]_| 


43.0m 


c* 


17.9± 4.0 


202 rpj 


12.23d 


c 


18. 9± 1.2 
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Continuation of Tabic 1 



Product 


Tl/2 


Type 


Yield (mb) 




72.912h 


c 


43.7± 2.9 




26.1h 


c 


40.6± 2.6 


200'j'2 


26.1h 


i(m+g) 


22.7± 1.5 




7.42h 


c 


38.5± 5.2 


198mlrpj 


1.87h 


i(ml+m2) 


17.6± 3.6 




5.30h 


c 


35.9± 5.0 


196mrj-ij 


84.6ni 


i(m) 


34.8± 4.4 


203 Hg 


46.612d 


c 


4.03± 0.27 


197m j_jg 


23.8h 


i(m) 


10.7± 0.7 


195m jjg 


41.6h 


i(m) 


13.6± 2.0 


193m jjg 


ll.Sh 


i(m) 


18.9± 2.5 


192 jjg 


4.85h 


c 


35.2± 2.8 


198m^y 


54.48h 


i(m) 


l.Olib 0.14 


198 Au 


64.684h 


i(m+g) 


2.11± 0.22 


198 Au 


64.684h 


i 


1.09± 0.30 


196Au 


6.183d 


i(ml+m2+g) 


4.13± 0.35 


195 Au 


186.098d 


c 


48.7± 5.5 


194Au 


38.020h 


i(ml+m2+g) 


7.06± 0.75 




4.94h 


c 


46.9± 6.6 


192Au 


4.94h 


i(ml+m2+g) 


11.6± 1.7 


191pt 


69.6h 


c 


40.1± 4.4 


189pt 


10.87h 


c 


46.8± 4.8 


188pt 


10.2d 


c 


40.5± 2.9 


190jj. 


11.78d 


c 


0.69± 0.06 


188 jj. 


41.5h 


c 


43.2± 3.2 


188jj. 


41. 5h 


i(m+g) 


2.93± 0.69 


186 jj. 


16.64h 


c* 


20.8± 1.9 


185 jj. 


14.4h 


c* 


34.8± 2.3 


1841^ 


3.09h 


c* 


39.5± 3.0 


1850s 


93.6d 


c 


41.8± 2.8 


183mOs 


9.9h 


i(m) 


23.2± 1.5 




22. Ih 


c 


42.0± 2.8 


183Re 


70.0d 


c 


41. 7± 2.9 




12.7h 


c 


45.2± 3.7 




19.9h 


c 


43.1± 5.9 




19.7m 


c* 


47.8± 4.2 


177W 


2.25h 


c 


30.1± 3.5 


176^ 


2.30h 


c 


30.8± 4.3 


i76Ta 


8.09h 


c 


34.5± 3.6 




3.14h 


c 


31.0± 3.9 


i72Ta 


36.8m 


c* 


17.3± 2.3 


175jjf 


70.0d 


c 


31.3± 2.3 


173Hf 


23.6h 


c 


28.4± 2.6 


172jjf 


683.017d 


c 


24.1± 1.6 


17lHf 


12.1h 


c 


18.2± 2.8 


170jjf 


16.01h 


c 


22.1± 6.8 


"2Lu 


6.7d 


c 


23.9± 1.7 


172Lu 


6.7d 


i(m+g) 


0.19± 0.05 


"iLu 


8.24d 


c 


26.1± 1.8 


i™Lu 


48.288h 


c 


21.7± 2.9 
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Continuation of Tabic 1 



Product 


Tl/2 


Type 


Yield (mb) 


i^^Lu 


34.06h 


c 


18.6± 1.2 


169Yb 


32.026d 


c 


20. 9± 1.5 


166Yb 


56.7h 


c 


16. 1± 1.1 




9.25d 


c 


19.4± 4.0 


165 rjj^ 


30.06h 


c 


14.4± 1.4 




28.58h 


c 


8.8± 0.6 




8.14h 


c 


5.73± 0.45 


155 


9.90h 


c* 


3.66± 0.27 


155Tb 


5.32d 


c 


4.16± 0.39 




56.16h 


0* 


2.52± 0.25 


152Tb 


17.50h 


c* 


2.10± 0.17 


153 Gd 


241. 6d 


c 


2.60± 0.23 


"9Gd 


9.28d 


c 


2.24± 0.18 


146 Gd 


48.27d 


c 


1.26± 0.09 


"^Eu 


24.0d 


c 


0.98± 0.31 


i^^Eu 


4.59d 


c 


1.63± 0.11 


146Eu 


4.59d 


i 


0.37± 0.05 




265.0d 


c 


1.02± 0.13 


i39Ce 


137.64d 


c 


0.83± 0.06 


121mrp^ 


154.0d 


i(m) 


0.44± 0.04 


i2iTe 


16.78d 


c 


l.llit 0.11 


119m "Jig 


4.7d 


i(ni) 


0.40± 0.04 


120msb 


5.76d 


i(m) 


0.54± 0.05 


114mjjj 


49.51d 


i(ml+m2) 


0.95± 0.19 


110m 


249.79d 


i(ni) 


1.12± 0.09 


106m^g 


8.28d 


i(m) 


0.89± 0.08 


105 Ag 


41.29d 


c 


0.65± 0.12 


105 Rh 


35.36h 


c 


4.63± 0.54 


loi'^Rh 


4.34d 


i(m) 


1.29± 0.16 


103Ru 


39.26d 


c 


3.84± 0.26 


96Tc 


4.28d 


i(m+g) 


1.20± 0.09 


95Tc 


20.0h 


c 


1.38± 0.13 


96 Nb 


23.35h 


i 


2.31± 0.19 


95Nb 


34.975d 


c 


5.41± 0.34 


95Nb 


34.975d 


i(m+g) 


3.03± 0.20 


95Zr 


64.02d 


c 


2.34± 0.15 


892r 


78.41h 


c 


2.30± 0.16 


882r 


83.4d 


c 


0.76± 0.08 


90m Y 


3.19h 


i(ni) 


4.82± 0.39 


88y 


106.65d 


c 


4.03± 0.27 


88y 


106.650d 


i(m+g) 


3.41± 0.25 


87y 


79.8h 


c* 


2.94± 0.23 


85sr 


64.84d 


c 


2.76 ± 0.22 


86Rb 


18.631d 


i(m+g) 


5.48 ± 0.66 


83Rb 


86.2d 


c 


3.46 ± 0.28 


S2™Rb 


6.472h 


i(m) 


2.73 ± 0.30 




35.3h 


i(m+g) 


2.17 ± 0.14 


75Se 


119.77d 


c 


1.34 ± 0.09 




17.77d 


i 


1.86 ± 0.18 


59pe 


44.503d 


c 


0.91 ± 0.08 


65Zn 


244. 26d 


c 


0.79 ± 0.19 




83.81d 




0.3o ± ().0() 
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Table 2: Experimental product nuclide yields in 2.6 GeV proton-irradiated ""*W 



Product 


Tl/2 


Type 


Yield (mb) 


ivvw 


2.25h 


c 


13.9 ± 1.9 


176^ 


2.30h 


c 


9.9 ± 2.9 




8.7h 


c 


4.44 ± 0.43 


i83Ta 


5. Id 


c 


10.5 ± 1.0 




114.43d 


c 


12.9 ± 1.3 


178mxa 


2.36h 


i(m) 


8.1 ± 1.3 


i^^Ta 


8.09h 


c 


29.3 ± 3.3 


i^^Ta 


10.5h 


c 


26.0 ± 2.8 


i^^Ta 


63m 


c 


25.8 ± 2.8 


ISljjf 


42.39 


c 


1.26 ± 0.12 


173Hf 


23.6h 


c 


29.9 ± 2.5 


17lHf 


12. Ih 


c 


19.6 ± 2.4 


170Hf 


le.oih 


c 


19.6 ± 4.0 




6.7d 


i(m-l-g) 


4.32 ± 0.56 




8.24d 


c 


30.1 ± 2.46 




8.24d 


i(m+g) 


10.8 ± 2.0 


i™Lu 


48.288h 


c 


24.8 ± 2.2 




34.06h 


c 


22.2 ± 1.8 




51.5m 


c 


23.4 ± 2.4 


167Yb 


17.5m 


c 


24.9 ± 2.8 


166 Yb 


56.7h 


c 


24.6 ± 2.1 


Tm 


7.7h 


c 


27.2 ± 2.3 


i66Tni 


7.7h 


i 


2.36 ± 0.46 


Tm 


30.06h 


c 


27.1 ± 2.4 


163 Tm 


l.Slh 


c 


26.3 ± 3.3 


i6iTm 


33m 


c 


21.0 ± 2.5 


leiEr 


3.21h 


c 


24.3 ± 2.5 


i60Er 


28.58h 


c 


23.9 ± 2.2 


i57Er 


25m 


c 


26.4 ± 5.9 


i56Er 


19.5m 


c 


15.9 ± 2.4 


160to1jJq 


5.02h 


i(ml+m2) 


24.9 ± 2.3 




12.6m 


c 


26.1 ± 6.8 


156Ho 


56m 


c 


19.6 ± 1.8 




8.14h 


c 


24.2 ± 2.2 


155j)y 


9.90h 


c 


22.1 ± 1.9 


153Dy 


6.4h 


c 


14.0 ± 1.9 




2.38h 


c 


15.6 ± 1.3 


155rpb 


5.32d 


c 


22.7 ± 1.9 


153Tb 


56.16h 


c 


18.9 ± 1.7 


152Tb 


17.50h 


c 


16.2 ± 1.3 




17.609h 


c 


16.7 ± 1.4 


149Tb 


4.118h 


c 


6.85 ± 0.62 


147Tb 


1.70h 


c 


2.15 ± 0.34 


i5iGd 


124.0d 


c 


19.0 ± 2.2 


149Gd 


9.28d 


c 


20.4 ± 1.7 


"7Gd 


38.1h 


c 


18.6 ± 1.6 


146Gd 


48.27d 


c 


19.4 ± 1.6 


"5Gd 


23.0m 


c 


12.9 ± 1.4 


149Eu 


93.1d 


c 


26.7 ± 3.4 



8 



Continuation of Tabic 2 



Product 


Tl/2 


Type 


Yield (mb) 


i^^Eu 


24.0d 


c 


22.4 ± 2.0 


146Eu 


4.59d 


c 


23.0 ± 1.9 


146Eu 


4.59d 


i 


3.62 ± 0.31 




5.93d 


c 


17.8 ± 1.6 




5.5h 


c 


2.87 ± 0.43 


i39Ce 


137.64d 


c 


19.8 ± 1.6 


i35Ce 


17.7h 


c 


17.8 ± 1.5 


CC 


3.51h 


c 


16.3 ± 2.7 




4.8h 


c 


14.5 ± 1.6 




11.50d 


c 


16.2 ± 1.3 


i26Ba 


100m 


c 


7.9 ± 1.1 




32.06h 


c 


18.7 ± 1.6 




36.4d 


c 


15.4 ± 1.3 


125XC 


16. 9h 


c 


14.2 ± 1.2 




2.08h 


c 


15.6 ± 1.3 




20.1h 


c 


11.7 ± 1.0 


i2iTc 


16.78d 


c 


10.7 ± 1.1 


ii9Te 


16.03h 


c 


9.17 ± 0.74 


119mr|-ig 


4.7d 


i(m) 


1.97 ± 0.17 




62m 


c 


8.81 ± 0.77 




5.0h 


i(m) 


1.08 ± 0.22 


"^Sb 


32.1m 


* 


9.85 ± 0.88 


ii3Sn 


115. 09d 




7.55 ± 0.67 


i"ln 


2.8049d 




7.44 ± 0.74 


llOm jjj 


4.9h 


i(m) 


3.29 ± 0.29 


109jn 


4.2h 


c 


5.12 ± 0.43 


106m^g 


8.28d 


i(m) 


1.70 ± 0.16 


105 Ag 


41.29d 


c 


5.33 ± 0.69 


lOOpd 


87.12h 


c 


1.24 ± 0.27 




20. 8h 


c 


3.97 ± 0.44 




20.8h 


i 


2.68 ± 0.28 


99m J^l^ 


4.7h 


c 


2.41 ± 0.28 


9^Ru 


69. 6h 


c 


3.13 ± 0.28 


96Tc 


4.28d 


i(m+g) 


1.73 ± 0.20 


93™Mo 


6.85h 


i(m) 


1.61 ± 0.13 


90Nb 


14.6h 


c 


2.58 ± 0.22 


892r 


78.41h 


c 


3.46 ± 0.28 


88 2r 


83.4d 


c 


2.56 ± 0.27 


88y 


106. 65d 


c 


3.49 ± 0.34 


SSy 


106.65d 


i(m+g) 


1.56 ± 0.22 


87y 


79. 8h 


c 


4.13 ± 0.34 


83Sr 


32.41h 


c 


1.96 ± 0.93 


84Rb 


32.77d 


i (m+g) 


1.31 ± 0.14 


83Rb 


86. 2d 


c 


3.34 ± 0.58 


82™Rb 


6.47211 


i(m) 


1.89 ± 0.17 




74.4m 


c 


1.71 ± 0.18 




119. 77d 


c 


2.38 ± 0.22 




7.15h 


c 


1.03 ± 0.11 




17.77d 


c 


1.38 ± 0.16 


Zn 


13.76h 


i(m) 


0.42 ± 0.038 


54 Mn 


312. 12d 


i 


2.51 ± 0.42 


51 Cr 


27.704(1 




4.5 ± 1.4 
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Continuation of Table 2 



Product 


Tl/2 


Type 


Yield (mb) 


4«v 


15.973d 


c 


0.557 ± 0.062 




43.67h 


i 


0.668 ± 0.091 


43k 


22.3h 


c 


0.681 ± 0.084 


28Mg 


20.91h 


c 


0.91 ± 0.089 


24Na 


14.95911 


c 


4.09 ± 0.34 


^Be 


53.29d 


i 


8.7 ± 1.0 



Simulation of experimental results 



Simulation techniques are of essential importance when forming the set of nuclear constants 
to be used in designing the ADS facilities, because they are universal and save much time and 
labour. At the same time, the present-day accuracy and reliability of the simulated results are 
inferior to experiment. Besides, the simulation codes are of different abilities to work when used 
to study the reactions that are of practical importance. 

Considering the above, the present work is primarily aimed at verifying the simulation 
codes used most extensively for the above purpose with a view to not only estimating their 
ability to work when applied to the issues discussed here, but also opening up ways to improve 
them. 

The following eight simulation codes were examined to meet these requirements: 

• the CEM95 cascade-exciton code [0], 

• the CASCADE cascade-evaporation-fission-transport code |^], 

• the INUCL cascade-preequilibrium-evaporation-fission code [0], 



the LAHET (ISABEL and Bertini options) cascade-evaporation-fission code [10 
the YIELDX semi-phenomenological code |11], 



the CASCADE/INPE cascade-preequilibrium-evaporation-fission-transport code |12 
the CEM2k cascade-exciton code p6|, a last modification of the CEM95 code, 



Contrary to the simulation results, the experimental data include not only the independent, 
but also (and mainly) cumulative and the supra cumulative yields of residual product nuclei. To 
get a correct comparison between the experimental and simulation results, theoretical cumulative 
yields must be calculated on the basis of the simulated independent yields. 

Since any branched isobaric chain can be presented to be a superposition of a few linear 
chains, the simulated cumulative and supra cumulative yields of a n-th nuclide can be calculated 

as 

n—l n—1 
„cum „ind < \ ^ And I I 



(7) 



i=l 



n-2 



aum ind , 

I7„ = cr„ + 



-I'n-l X 



a. 



^ind 
n — l 



(8) 
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Table 3: The yields (mb) of some products in the 1 GeV proton-irradiated ^ospj-, inferred from 
measurements at different laboratories; the ZSR and GSI data are taken from [6], the ITEP 
data are our present results 



Product nuclide 


ZSR Hannover 


ITEP 


GSI Darmstadt 


200 rp J 


22.3± 6.1 


22. 7± 1.5 


17.0± 0.4(1.6) 


iy«Au 


3.88± 0.47 


4.13± 0.35 


4.0± 0.1(0.4) 


ly^Au 


6.85± 0.92 


7.06± 0.75 


6.3± 0.2(0.6) 


l4«Eu 


0.104± 0.04 




0.075± 0.005(0.010) 


144pj^ 


0.068± 0.013 




0.036± 0.003(0.006) 



The branching ratios of the decay chains were retrieved from [13|. To get a correct comparison 



between results by different codes, the calculations were renormalized to unified cross sections 



for proton- nucleus inelastic interactions from [14|. 



If an experiment-simulation difference of not above 30% (0.77< Ucaic/(^exp <1.3) is taken 
to be the coincidence criterion |15|, the simulation accuracy can be presented to be the ratio of 
the number of such coincidences to the number of the comparison events. The 30% level meets 
the accuracy requirements of the cross sections for nuclide production to be used in designing 



the ADS plants, according to [15|. The mean simulated-to-experimental data ratio can be used 



as another coincidence criterion: 

< ^ > = 10\/< {o'cal,i/o-exp,i)'^ >^ (9) 

with its standard deviation 

S{< F>) = < (log {acai,i/<^exp,i) - log(< F >))2 >, (10) 

where <> designates averaging over all number of the experimental and simulated results 
used in a comparison. 

The mean ratio < F > together with its standard deviation S{< F >) defines the interval 
[< F >: S{< F >) , < F > xS{< F >)] that covers about 2/3 of the simulation-to-experiment 
ratios. 

The two criteria are considered sufficient to derive conclusions about the predictive power 
of a given code. The default options were committed to practical usage of the simulation codes. 



Comparison of data with simulation results 

The results obtained with the above-mentioned codes are presented in: 

• Figs, [s] and |6|, that show results of a detailed comparison between the simulated and 
experimental radioactive product yields; 

• Figs. and that show the simulated mass distributions of reaction products together 
with the measured cumulative (and supra cumulative) yields of the products that are at an 
immediate proximity to the stable isobar of a given mass (the sum of such yields from either 
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Table 4: Comparison statistics for ^^Pb and "'^*W 





Pb, E 


p = 1.0 GeV 


W, Ep = 2.6 GeV 


Code 


Nt = 116, Ng = 


= 95 


Nr = 107, Ng = 


= 93 




NC1.3/NC2.0/ 


< F > 


S{<F>) 


Nc,.3/Nc2.o/ 
/N5 


< F > 


S{< F >) 


LAHET 


41/65/90 


2.06 


1.78 


13/51/90 


2.52 


1.83 


CEM95 








28/66/81 


2.51 


2.23 


CEM2k 


38/58/66 


1.62 


1.44 


17/60/84 


2.24 


1.73 


CASCADE 


33/60/86 


2.28 


1.90 


48/71/91 


2.24 


2.04 


CASCADE/INPE 


36/66/84 


1.84 


1.56 








INUCL 


29/54/90 


2.87 


2.16 


38/58/86 


3.78 


3.23 


YIELDX 


30/54/90 


2.87 


2.24 


25/60/93 


2.04 


1.58 



sides in case both left- and right-hand branches of the chain are present). Obviously, the 
displayed simulation results do not contradict the experimental data if calculated values 
run above the experimental data and follow a general trend of the latter. This is because 
the direct 7-spectroscopy method used here identifies only radioactive products, that, as 
a rule, represents a significant fraction of the total mass yield, but, should a stable isobar 
of the given mass be produced, the 7-spectroscopy data are never equal to the total mass 
yield; 

• In Fig. m, that shows the experimental and simulated independent yields of reaction 
products in the form of isotopic mass distributions for several elements. 

Table ^ presents the statistics of our comparison between the experimental and simulated 
reaction product yields in the thin 208p|-, samples irradiated by 1.0 GeV and 2.6 GeV 

protons, respectively. Namely, it shows the total number of measured yields, N^^; the number of 
the measured yields selected to compare with calculations, Nq; the number of the product nuclei 
whose yields were simulated by a particular code, N5; the number of the comparison events when 
the simulated data differ from the experimental results by not above 30%, Nci.s; the number 
of the comparison events when the simulated data differ from the experimental results by not 
more than a factor of 2.0, Ng2.o- 

Since about 30% of all measured secondary nuclei are not spallation reaction products, an 
important criterion of the codes is their ability to simulate the high-energy fission and fragmen- 
tation processes. Among the codes used here, LAHET, CASCADE, INUCL, CASCADE/INPE, 
and YIELDX simulate both spallation and fission. The CEM95 and CEM2k codes simulate 
spallation only, which is explicitly reflected in a smaller number of the products simulated (the 
parameter N^ in Table 4 and in the shapes of the simulation curves in Figs. 5-8. 

The following conclusions follow from our analysis of the experiment-to-simulation com- 
parison results presented in Table 4 and in Figs. |5[|^: 

1. Generally, all codes can quite adequately simulate the weak spallation reactions (the A>180 
products for ^'^'^Pb and the A>150 products for "'**W), with the simulation results differing 
from experimental data within a factor of 2. 
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2. In the deep spallation region (150<A<180 for ^ospb and 110<A<150 for "'^*W), the sim- 
ulation codes are of very different predictive powers, namely, 

• the LAHET (when not shown explicitly as "Bertini" , all results by LAHET are of the 
ISABEL option), CEM2k, CASCADE/INPE, and YIELDX predictions are actually 
the same as the experimental data; 

• the CASCADE code simulates the A>160 product yields adequately. Below A = 160, 
however, the simulated data get underestimated progressively (up to a factor of 5) 
compared with experiment (see Fig. 7); 

• the INUCL code underestimates the yields of all the products by a factor of 2-10 in 
all the above mass ranges (see Figs. 7 and 8). 

3. In the mass range characteristic of the fission products (50<A<150 for 208p|-, ^nd 30<A<110 
for ""*W), the INUCL code predictions are in the best agreement with experiment when 
describing the yields from ^'^^Pb. As a rule, the INUCL-simulated results differ from the 
data by not above a factor of 1.5. In the case of "'^*W, however, the prediction quality 
deteriorates substantially. The LAHET-simulated yields are underestimated by a factor 
of 1.5-10.0 for Pb (Figs. 5 and 7) in the whole fission product mass region and for A<60 
in the case of W (Figs. 6 and 8) but are oversetimated several times for fission fragments 
with A>60 from W. The YIELDX-simulated yields are either under- or over-estimated by 
a factor of up to 30 without showing any physical regularities. The CASCADE/INPE- 
simulated yields of the 130<A<150 reaction products are strongly underestimated (up to 
1-2 orders of magnitude), while the simulated 40<A<130 product yields agree with the 
data within a factor of 2, as a rule. Generally, all the codes exhibit the feature noted above 
for INUCL, namely, the yield prediction quality in the case of "''^W is much worse com- 
pared with ^''^Pb, probably, because the fission cross sections of high-excited compound 
nuclei with very low fissility are difficult to calculate. 



4. The last version of the improved cascade-exciton model code, CEM2k |16|, shows the best 



agreement with the 1 GeV Pb-data in the spallation region, especially for the isotopic mass 
distributions (Fig. 9). At 2.6 GeV (W-target), it overestimates the expected experimental 



fission cross section of about 41 mb |17] by a factor of 6. This overestimation of the fission 
cross section causes an underestimation of the yield of nuclei which are most likely to 
fission (with a very low fissility) at the evaporation stage of a reaction, after the cascade 
and preequilibrium stages , i.e., for 147 < A < 175 (see Fig. 6). Similar disagreement with 
the 2.6 GeV W-data one can see as well for LAHET and CEM95, that is also related with 
an overestimation of the fission cross section at 2.6 GeV (see Figs. 6 and 8). The code 
CEM2k is still under development, its problem with the overestimation of fission cross 
sections at energies above 1 GeV has yet to be solved, and it has to be complemented with 
a model of fission fragment production, to be able to describe as well fission products. 



Conclusion 



The trends shown by the advances in the nuclear transmutation of radioactive wastes 
and Spallation Neutron Source (SNS) facilities permit us to expect that the accumulation and 
analysis study of nuclear data for ADS facilities will have the same rise of academic interest and 
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Figure 5: Product comparison between the experimental (closed symbols) and simulated (open 
symbols) yields of radioactive reaction products from 208p|-, irradiated with 1 GeV protons. The 
cumulative yields are labeled with a "c" when the respective independent yields are also shown. 



■ EXPERIMENT O YIELDX<F>=2.04 




Figure 6: Product comparison between the experimental (closed symbols) and simulated (open 
symbols) yields of radioactive reaction products from '^"^W irradiated with 2.6 GeV protons. 
The cumulative yields are labeled with a "c" when the respective independent yields are also 
shown. 
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Product mass number 

Figure 7: Measured and calculated by the codes mass product yields from ^osp^ irradiated with 
1.0 GeV protons. For comparison, the GSI data from [6] are shown as well. 



practical commitments as in the nuclear reactor data during the last five decades. Therefore, 
the experimental data on the yields of the proton-induced reaction products as applied to the 
ADS and SNS main targets and structure materials are urgent to accumulate. It should be 
emphasized that the charge distributions in the isobaric decay chains are important to study as 
well. The data thus obtained would make it possible, first, to raise the information content of the 
comparisons between the experimental and simulated results and, second, to lift the uncertainties 
in experimental determination of the cumulative yields by establishing unambiguous relations 
between cj'^"™- and cj'^"™- for many of the reaction product masses. 

Regarding the codes benchmarked here, one may conclude that none of them agree well with 
the data in the whole mass region of product nuclides and all should be improved further. The 
new CEM2k code developed recently at Los Alamos [16| agrees with our data in the spallation 
region the best of the codes tested. But CEM2k has yet to be completed by a model of fission 
fragmentation, to become applicable in the fission-product region as well. 
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Figure 8 
2.6 GeV 



Measured and calculated via the codes mass product yields from ""*W irradiated with 
protons. 
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